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Employing  a  semi-classical  model,  numerical  solutions 

were  obtained  for  a  ten  level  model  of  methyl  fluoride. 

Laser  emission  from  methyl  fluoride  at  approximately 

16  microns  was  found  to  be  possible  by  optically 

pumping  with  a  carbon-dioxide  laser.   The  methyl  fluoride 

laser  was  found  to  have  gain  of  approximately  .51%  per 

cm.  and  energy  storage  of  approximately  31  mJL" l torr" l . 


CHAPTER  I  -  INTRODUCTION 

Laser  (La'zSr)  n.  [l(ight)  a (mplif ication  by) 
s(timulated)  e  (  mission  of  )  r  (  adiation)  ]  ;  a  device, 
containing  a  crystal,  gas,  or  other  suitable  substance,  in 
which  atoms,  when  stimulated  by  focused  light  waves, 
amplify  and  concentrate  these  waves,  then  emit  them  in  a 
narrow,  very  intense  beam,  optical  maser. x 

Thus  Webster  defines  the  laser.   Webster's  definition 
exemplifies  the  importance  and  rapid  growth  of  this  new 
light  source.   Prior  to  the  summer  of  1960  however   the 
word  laser  had  an  entirely  different  meaning;  "Laser;  the 
juice  of  the  laser  tree  ,   laserpitium  latifolium,  also 
called  silphium,  greatly  esteemed  by  the  ancients  as  an 
antispasmodic  ,  deobstrucent  and  diuretic."2   Although 
the  word  has  an  ancient  heritage,  today  it  brings  to  mind 
a  much  more  modern  meaning.   The  laser  is  an  example  of 
the  same  idea  ocurring  to  different  people  in  different 
parts  of  the  world  at  about  the  same  time.   The  first 
lasers  were  actually  masers  (microwave  amplifiers)  which 
were  jointly  developed  in  1951  by  C.  H.  Townes  of  Columbia 
University  and  1952  by  N.  G.  Basov  and  A.  M.  Prokhorov  of 
the  USSR.   For  their  work  on  coherent  microwave  amplifiers 
Townes,  Basov,  and  Prokhorov  were  jointly  awarded  the  Nobel 
prize  in  Physics  in  1964.   Following  the  development  of  the 
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maser,  A.  L.  Schawlow  of  Bell  Telephone  Laboratories  and 
Townes  outlined  the  theory  of  an  optical  maser  in  1958. 
After  publication  of  this  theory,  T.  H.  Maiman  in  1960  was 
the  first  to  obtain  actual  laser  action  in  the  visible 
spectrum.   Maiman  employed  a  crystal  of  ruby  as  a  lasing 
medium.   This  was  suprising  since  at  that  time  scientists 
were  working  on  primarily  gas  lasers,  and  felt  that 
probability  of  success  was  highest  for  gaseous  devices. 
Time  has  proven  their  premonitions  correct  since  most  high 
power  laser  development  in  recent  years  has  been  from 
gases.    Solid  state  lasers  have , however , been  found  to  be 
sources  of  tunable  radiation  important  to  ultra-high 
resolution  spectroscopy,  and  recently  glass  lasers  have 
been  used  in  controlled  fusion  studies.   Outside  of  these 
two  areas  new  high  power  laser  applications  are  being 
developed  primarily  around  gas  lasers. 

Prior  to  the  late  sixties  the  laser  was  largely  a 
scientific  toy.   It  was  at  that  time  that  molecular  gas 
lasers  were  developed.   These  lasers , particularly  the  C02, 
have  exhibited  several  important  qualities.   First  is 
their  efficiency  (30-40%) .   This  is  extremely  high  when 
compared  with  efficiencies  of  less  than  1%  for  most 
conventional  lasers.  Second,  and  perhaps  more  important, 
is  their  ability  to  be  scaled  up  to  industrial  size. 
This  gave  great  impetus  toward  developing  commercially 
available  high  power  C02  systems.   In  recent  years  C02 
lasers  have  been  moving  out  of  the  scientific  laboratory 


3 
and  into  industry  and  the  production  line.3'^'5 

Continuous  wave  (CW)  lasing  at  10.6  microns  from  C02 
was  first  reported  by  C.  N.  K.  Patel  in  1964. 6   Since  that 
date  the  C02  laser  has  undergone  rapid  improvement. 
Attention  is  called  to  the  review  of  C02  laser  development 
by  Robinson. 8 

Laser  action  from  C02  proceeds  via  population  inversion 
of  vibrational  levels.   In  electric  discharge  systems 
nitrogen  is  excited  by  electron  collisions,  which  in  turn 
excites  the  00 ° 1  antisymmetric  stretch  mode  (see  figure  1) 
of  C02.   Lasing  is  achieved  from  this  level  to  either  the 
10 °0 symmetric  stretch  at  10.6  microns  or  the  02 °0  second 
bending  mode  at  9.6  microns  (figure  2).9   In  order  to 
maintain  a  population  inversion  the  lower  (01*0)  level  must 
be  depleted  .  Normally  helium  is  added  to  collisionally 
depopulate  this  level.   In  addition  to  undergoing  ir  active 
vibrational  energy  transitions,  C02  may  change  its  rotational 
energy  by  one  quantum.   This  leads  to  a  vibration-rotation 
spectrum  of  many- frequency  laser  lines  in  each  the  10.6 
micron  band  (figure  3) .    Each  term  labeled  by  J  represents 
a  rotational  angular  momentum  state  of  C02  and  thus  a 
splitting  of  energy  levels.10  The  terms  R(J)  and  P (J) 
represent  a  transition  of  AJ  =  -1  and  +1  respectively, 
following  ir  selection  rules  for  C02 .   Figure  4  illustrates 
the  wavelength  difference  for  P,Q,  and  R  bands  (Q;AJ  =  0). 
Notice  that  odd  J  terms  are  missing.   This  is  due  to 
nuclear  spin  statistics  for  zero  nuclear  spin  molecules 
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Figure    1.    Energy   levels   of  C02 
and   N2    laser    levels. 
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Figure  2.  Laser  vibrations  of 
C0o. 
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Figure  3.  V-R  spectrum  of  CO- 
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Figure  3.  V-R  spectrum  of  CO  . 


Figure  4.  Diagram  of  V-R  selection  rules, 
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(for  further  discussion  see  Herzberg) . 

The  theory  of  laser  oscillation  is  easily  understood 
if  presented  within  the  framework  of  amplifier  dynamics. 
The  laser  as  a  device  consists  usually  of  a  pair  of 
parallel  mirrors  between  which  is  placed  the  material  which 
is  to  act  as  the  amplifying  medium  over  a  limited  frequency 
range  (figure  5) .   L  is  the  length  of  the  active  material 
and  the  optical  distance  will  be  denoted  L".   If  mirrors 
are  placed  on  the  ends  of  the  cavity  containing  the 
amplifying  material  L'=  nL  where  n  is  the  index  of 
refraction.   For  gas  lasers  n-1   however  the  mirrors  are 
usually  located  outside  the  laser  cavity  so  that  L">  L. 
At  best  one  of  the  reflectors  is  partially  transmitting  so 
that  its  reflection  coefficient  (reflectivity)  is  less 
than  one.   The  reflection  coefficient  r   is  defined  as  the 
fraction  of  light  intensity  reflected.   Thus  at  each 
reflection,  (1-r)  of  the  intensity  is  transmitted  as  output. 
If  the  reflectivities  of  the  two  mirrors  are  ri  and  r2 
respectively,  the  energy  of  the  wave  is  diminished  by  a 
factor  rxr2  in  one  transit.   If  the  fraction  of  intensity 
remaining  after  one  round  trip  through  the  laser  is  denoted 
by  e2^  ,  Y  is  seen  to  be  the  loss  coefficient  per  length. 
It  is  seen  that 

y    =   -1/2  Log  rir2 
if  all  other  losses  other  than  reflection  are  neglected. 
In  gas  lasers,  diffraction  losses  are  sometimes  quite 
significant,  but  they  will  be  neglected  here. 
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Figure  5.  Diagram  of  gas  laser. 
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Oscillation  may  be  sustained  if  the  amount  of  amplifi- 
cation per  passage  is  equal  to  or  greater  than  the  losses. 
If  the  amplification  per  unit  length  is  denoted  by  am, 
the  threshold  condition  is  seen  to  be 

(1)  amL  =  y. 

This  equation  also  defines  the  amplification  since  the  laser 
will  always  operate  at  threshold,  as  in  electrical  circuits. 
The  amplification  coefficient  (negative  of  the  absorption 
coefficient)  may  be  written  as 

(2)  a(v)  -  k(v)N 
where  v  =  frequency 

N  =  relative  population  difference  between  emitting 
and  absorbing  levels 

N  =  1/N0  {qjq2    Nj-Nj 
g1 ,    g2  are  the  quantum  degeneracies  of  levels  1  and  2 
respectively-   N   and  N2  are  the  populations  of  levels  1 
and  2  respectively.   It  can  be  shown  that 

(3)  k(v)  =  a?1  /8^n2v2  g2/g1  G(vv°)  =  kG(w°) 
where  A   =  Einstein  spontaneous  emmission  coefficient 

(to  be  discussed  later) 
A2~\l   =  t2  the  lifetime  of  atoms  in  level  2. 
H  =  index  of  refraction 
G(vv°)  =  lineshape  around  line  center  v°. 
Combining  equations  (1) ,  (2) ,  and  (3)  we  find 

(4)  g,/g2  N2-  Nj  =  y/LkG(0) 
where  we  assume  emmision  is  at  line  center. 
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The  time  for  a  photon  to  make  a  single  laser 
passage  is  t  =  I//C  so  that  P0  photons  travelling 
through  the  laser  after  m.  passages  will  be  reduced  to 
p  _  p  e"Yt/T 

thus  the  average  photon  lifetime  may  be  defined  as 

(5)     tp=  t/y  =  L'/cy. 
Combining  this  with  equation  (4) 


2,,  2 


i  r\  a ,  8im  v  g,  L   t2 

(6)      ftN2-  n,=  cwrltL  17  ' 

Equation  (6)  is  known  as  the  Schawlow-Townes  condition, 
and  gives  the  minimum  population  inversion  necessary 
for  lasing. 

Once  a  population  inversion  has  been  established  the 
principle  of  stimulated  emission  may  be  simply  stated  as: 
A  photon  will  interact  with  a  molecule  in  an  excited  state 
to  produce  two  coherent  photons  travelling  the  same  direct- 
ion as  the  first,  plus  a  molecule  in  the  ground  state. 
The  initiation  and  maintenance  of  lasing  action  is  most 
easily  explained  by  use  of  the  Einstein  coefficients. 
These  two  coefficients  are 
spontaneous  A12  [sec  1]    and 
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stimulated  B  l2  [cm3J  '  sec   ]  . 

An  excited  state  has  probability  Pmn  of  decaying  from 

state  m  to  state  n  defined  as 

pmn  =  Amn  +  uVBmn 
where   U   is  the  radiative  density  at  the  frequency 
satisfied  by  the  Bohr  condition.   Note  that  Amn  and  Behave 
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different  units.   Amn  is  a  measure  of  the  natural  life- 
time of  state  m  before  spontaneous  radiative  decay  to  state 
n,  analogous  to  noise  in  an  electrical  amplifier.   Just  as 
in  the  electrical  analogy,  the  spontaneously  emitted  photons 
start  the  oscillations  which  lead  to  a  radiation  density 
thus  stimulating  more  photons.   Once  the  oscillations  have 
reached  threshold,  the  stiumlated  term  will  continue  to 
act  as  an  amplifier  as  long  as  a  population  inversion  is 
maintained.   It  is  common  in  lasers  for  population  inversion 
to  be  established  and  oscillation  to  start,  where  the  rate 
of  depletion  of  the  upper  level  is  faster  than  it  can  be 
refilled.   In  this  case  gi/g?  Ni-  N2  <  0  within  a  short 
time,  and  laser  oscillations  cease.   This  type  of  laser  is 
called  a  pulsed  laser,  as  opposed  to  a  CW  laser.   C02  may 
be  used  as  either  and  may  be  constructed  in  several 
different  configurations.   The  earliest  and  simplest  was  by 
Patel  (see  figure  6) . n  It  is  interesting  to  note  that  the 
C02  was  not  placed  directly  in  the  discharge.   Nitrogen  is 
excited  in  the  discharge  and  being  homonuclear,  is  not  ir 
active.   C02  is  then  collisionally  excited  by  the 
nitrogen  to  achieve  a  population  inversion.   Today  commer- 
cial lasers  place  a  1-1-8  mixture  of  N2-  C02-  He  in  the 
discharge  to  obtain  pulsed  or  CW  lasing. 

An  important  advancement  in  C02  laser  technology  was 
the  development  of  the  TEA  (Transverse  Electric  Atmospheric) 
laser.12  As  the  acronym  suggests  this  laser  operates  at 
atmospheric  pressures  thus  eliminating  the  need  for  vacuum 
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Figure  6.  Early  CO2  laser  by  Patel, 


15 

systems  and  increasing  power  output.   The  pulse  length  is 
largely  determined  by  the  relaxation  times  of  the  gas 
molecules.   As  these  become  shorter  at  higher  pressures, 
the  pulse  length  decreases  proportionally.   Typical 
pulse  length  of  a  TEA  laser  is  a  few  nanoseconds 
(1  nsec  =  10~9sec).   It  is  difficult  to  achieve  uniform 
excitation  of  C02  at  atmospheric  pressure,  since  at  about 
200  Torr  the  normal  glow  discharge  constricts  to  an  arc, 
heating  the  gas  and  destroying  laser  action.   Several 
techniques  have  been  developed  to  overcome  this. 
Beaulieu   originally  used  a  very  short  17  KV  pulse  to 
excite  a  transverse  discharge  between  a  row  of  pin 
cathodes  and  a  round  bar  anode.12  Today  numerous   TEA 
lasers  use  a  helix  geometry  transverse  discharge.   The 
helix  is  designed  to  improve  the  uniformity  of  the  discharge. 

The  radial  radiation  distribution  in  a  cylindrical 
laser  is  very  close  to  the  eigenf unctions  of  a  cylinder 
(see  figure  7) .   Because  of  diffraction  losses,  the  TEMQ0 
mode  propagates  best,  and  to  prevent  mode  switching  and 
competition  most  lasers  are  operated  by  placing  a  stop  with- 
in the  optical  cavity  so  that  the  laser  operates  in  TEM00  . 

It  was  indicated  earlier  that  gas  lasers  radiate  on 
many  vibrational-rotational  transitions.   Because  cavity 
quality  factor  is  narrow  compared  with  the  Doppler  or 
pressure  broadened  half  width,  the  laser  has  longitudinal 
modes  within  the  natural  half  width.13  The  condition  for  re- 
inforcement is  that  nX=2Ln  or  v/c=n/2Ln ,  defining  the 
optical  distance 


TEM 
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Figure  7. Radial  distribution  of  radiation  for 
cavity  with  circular  mirrors.  TEM  (Transverse 
Electric  and  Magnetic)  modes  are  indicated  by 
arrows  indicating  magnitude  of  Electric  field. 
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1/  =  /0  n(z)  dz 
the  condition  for  reinforcement   becomes 

v/c  =  n/2L" . 

Thus  the  spacing  between  consecutive  axial  modes  is 

v  .   -  v   =  Av  =  c/2L'. 
vn+i     n 

Axial  mode  spacing  is  compared  with  natural  line  width  in 
figure  8,  notice   that  the  cavity  spacing  is  narrow  with 
respect  to  the  natural  line.   Only  those  cavity  modes  above 
threshold  will  lase,  and  in  most  gas  lasers  several  axial 
modes  will  operate  simultaneously  giving  multimode  operation. 
There  may  possibly  be  from  several  to  several  thousand  modes 
excited  during  laser  operation. 

It  is  of  interest  to  consider  the  type  of  output  that 
would  result  if  a  large  number  of  modes,  equidistant  in 
frequency,  could  be  excited  with  approximately  equal 
amplitude,  and,  if  these  modes  could  somehow  be  locked 
together  in  a  constant  phase  relationship.   Each  oscillating 
mode  gives  rise  to  a  wave  with  an  amplitude  described  by 
the  real  part  of  eia3(t_x/c)  where  u  =  w0  +  kAto  and  Aw  =  2ttAv 
(Av  =  half  width  of  mode) .   The  addition  of  2n  +  1  of  such 
waves  with  frequency  centered  around  w„  all  having  equal 
amplitude  leads  to  the  expression 

ln  ei(u)o  +  kAo))t 

k=-n 


which  is  equal  to e1Wot  times  the  amplitude  function 

F(t)  =  sin(n+^)  Acot.   The  result  of  this  synthesis  is  an 
sin  7Acot 

amplitude-modulated  wave  of  frequency  u>0  ,  whose  intensity 
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Figure  8. Axial  mode  spacing. 
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varies  at  the  rate  F(t)2.   Figure  9  illustrates  a  typical 
example.   The  peak  intensity  reaches  (2n+l) 2  and  the  first 
zero  of  F(t)  is  separated  from  the  peak  by  a  time  interval 
t  =  T/2n+l  where  T  is  the  period  equal  to  2tt/Aco.   For  a 
laser  of  optical  length  L'  =  60cm,  T  =  4xl0"9sec.   If  100 
adjacent  axial  modes  are  locked  together  the  resulting  peak 
intensity  is  10**  times  that  of  the  individual  modes,  that 
is,  100  times  the  sum  of  all  intensities.   The  peaks 
would  repeat  at  the  rate  off1  =  250MHz  and  peak  pulses 
would  have  a  half  width  of  4x10" n  sec. 

Experimentally,  phase-locking  of  the  axial  modes  was 
first  realized  by  means  of  an  acoustic  modulator  incorporated 
in  a  He  -  Ne  laser  .   Further  experimentation  with  mode- 
locking  revealed  that  it  could  be  accomplished  without  the 
use  of  an  externally  driven  modulator.   Incorporation  of  a 
suitable  bleachable  absorber  cell  produces  self-locking  of 
the  longitudinal  modes.   Figure  10  illustrates  schematically 
the  use  of  a  bleachable  absorbing  cell  as  a  mode  locker. 
The  complete  description  of  this  process  requires  a  detailed 
mathematical  analysis  although  several  things  can  be  said 
qualitatively.11*  The  modes  are  all  locked  in  phase  when  the 
dye  cell  is  short  and  is  located  near  one  of  the  mirrors. 
Mode  locking  is  accomplished  by  the  nonlinear  interaction 
of  radiation  of  differing  frequencies  within  the  dye  cell. 
The  nonlinearity  of  the  interaction  is  the  consequence  of 
the  fact  that  the  dye  cell  is  operated  in  a  region  of 
intensity  where  the  transitions  are  nearly  saturated. 
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Figure  9.  Example  of  mode  locked  pulses. 
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Figure    10.    Mode    locked    laser. 


22 
The  need  for  saturable  absorbers  was  one  of  the  early 
motivations  for  finding  gases  which  efficiently  absorbed 
laser  radiation.   In  the  case  of  C02  ,  several  gases  were 
found  such  as  SF6 ,  NH3,  PF5,  and  CH3F.     It  was  soon 
discovered  however  that  laser  absorption  provided  a  means 
of  preparing  a  molecule  in  an  excited  state  in  order  to 
study  energy  transfer  processes  and  chemical  reactions.   The 
criterion  for  good  absorption  of  laser  radiation  is  that  the 
gas  must  have  an  energy  level  resonant  with  the  laser 
emission.   For  infrared  lasers,  this  corresponds  to 
molecular  vibrational  frequencies,  and  it  was  soon  dis- 
covered that  the  laser  could  play  an  important  role  in  the 
study  of  a  process  mediated  through  vibrations  such  as 
chemical  reaction.16-23   Lasers  have  become  an  important  tool 
in  fluorescence  studies,  ultra  high  resolution  spectroscopy, 
and  molecular  energy  transfer,  as  well  as  developing  passive 
mode  lockers  for  other  lasers.2"  Figure  11  shows  a  typical 
experimental  set-up  for  a  quasi-  CW  fluorescence  experiment. 
Infrared  fluorescence  may  be  monitored  from  a  gas  which 
absorbs  C02  radiation  directly,  or  by  the  addition  of  a 
low  concentration  of  sensitizer  such  as  SF6 ,  the  fluorescence 
from  the  collisionally  excited  molecule  may  be  observed. 
Laser  induced  fluorescence  spectroscopy  is  also  a  valuable 
technique  for  following  the  course  of  laser  induced  chemical 
reactions.25'26  Molecules  can  be  excited  by  the  absorption 
of  C02  laser  radiation,  and  the  kinetics  of  the  subsequent 
reactions  can  be  monitored  by  ir  fluorescence  spectroscopy. 
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Figure  11.  Experimental  set-up  of  laser  fluorescence 

study. 
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In  general  the  reactions  are  monitored  by  observing  the 
fluorescence  from  a  given  vibrational  level  as  a  function 
of  time.   Detailed  energy  transfer  and  evidence  of  non- 
Boltzmann  chemical  reactions  have  been  observed  in  CH3F, 

CH3CI'  and  CH2C12.  ' 

If, as  postulated,  chemical  processes  can  be  affected 

by  selectively  exciting  vibrational  levels  in  one  of  the 
reactants  producing  non-equilibrium  distributions, 
preferential  reaction  channels  or  increased  reaction 
rates  may  be  obtained.   Possible  applications  of  this 
technique  include  accelerated  catalysis,  efficient 
fractionation  of  hydrocarbons,  and  isotope  separation. 
To  date  the  possibility  of  stimulating  chemical  reactions 
with  lasers  has  been  shown  to  be  possible;  however  as 
yet  there  has  been  no  large  scale  applications  of  non- 
Boltzmann  chemical  reactions.26 

The  possibility  of  separating  isotopes  with  lasers 
has  received  much  attention  in  recent  years.   The  central 
idea  is  selective  excitation  of  the  isotopic  species  by 
a  narrow  band  ir  laser.   This  occurs  because  of  the 
isotope  shift, which  is  in  the  infrared  portion  of  the 
spectrum,  so  that  the  laser  leaves  only  a  single  isotope 
vibrationally  excited.   Photoionization  may  proceed 
from  the  excited  molecule  by  a  visible  laser  which 
would  not  contain  enough  energy  to  ionize  the  non- 
excited  molecule.   The  ionized  molecule  may  then  be 
removed  electrostatically  or  chemically.23 
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One  major  problem  that  must  be  overcome  is  the 
exchange  of  charge  or  vibrational  energy  between  isotope 
and  non-isotope  before  the  separation  takes  place.   This 
would  mix  the  molecules  enough  so  that  no  substantial 
increase  in  isotope  concentration  would  result.   Lasers 
may  be  the  solution  since  a  laser  induced  reaction  might 
proceed  rapidly. 

Although  the  idea  of  laser  isotope  separation 
appears  sound,  in  practice,  few  isotopes  have  been 
separated  by  either  the  photoionization ,  photon  recoil, 
or  induced  chemical  reaction  method.   One  successful 
isotope  separation  has  been  with  the  isotope  of  clorine. 
The  reason  for  this  success  is  the  fact  that  BC13  has 
an  energy  level  resonant  with  C02  laser  radiation. 
Not  surprisingly, boron  isotopes  have  also  been  separated 
by  this  process.   There  are  however  many  isotopes  which 
do  not  absorb  resonantly  with  conventional  lasers.  It 
is  the  separation  of  these  isotopes,  which  include 
uranium, that  have  stimulated  interest  in  developing 
middle  infrared  lasers  of  high  (several  percent) 
efficiency. 

The  most  important  isotope  separation  problem 
today  is  the  U235  -  U 238  pair  for  obvious  reasons.   A 
glance  at  the  spectrum  of  UF6  reveals  that  the  absorp- 
tion at  16  microns  is  by  far  more  intense  than  other 
wavelengths  and  the  isotope  shift  is  largest  there  also. 
A  16  micron  gas  laser  does  not  currently  exist  and  its 
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development  has  stimulated  interest  in  several  laboratories. 

A  new  development  that  appears  promising  in  genera- 
tion of  ir  lasers  is  optical  pumping  with  a  conventional 
laser  source.   This  eliminates  the  need  for  operation  in 
an  electric  discharge  or  flowing  gas  chemical  system 
which  usually  requires  a  homonuclear  or  metastable 
collision  partner  such  as  N2  and  H2.   The  fundamental 
idea  is  absorption  of  laser  radiation  of  one  wavelength 
and  regeneration  at  another.   The  quantum  efficiency  of 
such  a  process  would  be  the  ratio  of  the  two  wave- 
lengths or  about  62%  for  C02  absorption  and  16  micron 
emission,  although  quantum  efficiency  is  seldom  approached. 
C02  lasers  commonly  operate  with  20%  efficiency  which 
would  yield  a  16  micron  laser  with  1-2%  efficiency. 
Although  middle  infrared  lasers  have  not  yet  been 
generated  in  this  manner,  nearly  200  different  wave- 
lengths lasers  in  the  far  infrared  have   and  it  is 
hoped  that  this  method  will  eventually  provide  lasers 
at  nearly  any  desirable  wavelength.  3(h  32  The  application 
of  optical  pumping  for  generation  of  middle  infrared 
lasers  to  be  used  in  the  study  of  isotope  separation, 
laser  induced  chemistry,  and  ultra  high  resolution 
spectroscopy  should  deserve  more  and  more  attention  in 

the  future. 

The  author  believes  an  important  problem  which 
may  be  solved  by  optical  pumping  is  the  development 
of  a  laser  at  approximately  16  microns  for  the  separation 
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of  uranium  isotopes.   By  using  a  conventional  laser 
source  such  as  C02 ,  a  molecule  may  be  excited  and  made 
to  regenerate  at  around  16  microns.   Figure  12  shows  a 
diagramatic  representation  of  the  energy  levels  required 
for  successful  absorption  of  C02  radiation  and  regenera- 
tion at  16  microns. 

The  purpose  of  this  paper  is  to  demonstrate  the 
feasibility  of  generation  of  middle  infrared  lasers 
by  optical  pumping.   The  particular  wavelengths  used 
will  be  absorption  of  C02  radiation  and  emission  at 
16  microns,  although  the  ideas  should  be  equally 
applicable  to  other  needs. 
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Figure  12.  Expected  energy  levels  leading 
to  a  16micron  laser. 


CHAPTER  II  -  ENERGY  TRANSFER  IN  MOLECULES 

Molecular  energy  transfer  is  important  to  both  laser 
absorption  and  laser  emission  processes.   For  radiation 
in  the  infrared,  energy  transfer  is  largely  in  vibra- 
tional and  rotational  modes.  Schrodinger* s  equation  is 
written  by  assuming  the  Born-Oppenheimer  approximation 
(4)     ¥(r,R)  =  ye(r)Xn(R) 

where  ¥  (r)  is  the  electon  wavefunction   and  xn  ^s  tnat 
for  the  nucleus.   x   includes  all  nuclear  terms  such  as 
vibrational  and  rotational  motion.   The  energy  is  the 
sum  of  the  individual  contributions  viz, 
Etot=  Etrans  +  Erot  +  Evib  +  Eel  +  Espin-orbit 

Normally  the  contributions  from  infrared  processes  are 
those  of  vibrational  and  rotational  motions.   The 
electronic  configuration  is  usually  considered  constant 
for  infrared  processes.   Eigenf unctions  for  vibrational 
motion  are  usually  found  by  first  assuming  a  harmonic 
potential  and  adding  anharmonic  perturbations.   The 
eigenfunctions  for  a  harmonic  potential  are  harmonic 
oscillator  wavefunctions  which  lead  to  an  evenly  spaced 
spectrum  with  selection  rules  Av  =  ±1  (see  figure  13) . 
The  addition  of  anharmonic  terms  perturbs  the  spectrum 
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Figure    13.    Spectrum  of  harmonic   oscillator, 
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qualitatively  as  shown  by  the  dashed  lines. 

As  indicated  earlier,  rotational  energy  splits  each 
vibrational  energy  level  into  many  sublevels,  the 
eigenfunctions  of  which  are  found  by  assuming  a  rigid 
rotor  Hamiltonian  and  are  proportional  to  spherical 
harmonic  functions.   This  leads  to  selection  rules 
A  J  =  0,  ±1,  for  a  totally  symmetric  electronic  wave- 
function.33   Figure  14  shows  the  energy  levels  of  a 
rigid  rotor  and  indicates  dipole-allowed  transitions. 

Emission  or  absorption  of  light  is  usually  explained 
within  the  framework  of  time  dependent  perturbation 
theory.31*   After  separating  the  unperturbed  Hamiltonian 
into  vibrational,  rotational , and  electronic  motion, 
solutions  are  found  in  the  absence  of  any  electromagnetic 
radiation.   The  oscillating  field  is  then  treated  as  a 
perturbation  and  first  order  quantities  such  as  probabil- 
ities of  transition  from  one  stationary  state  to  another, 
and  energy  corrections  are  calculated.   If  one  writes 
the  perturbing  hamiltonian  as 
H'(Rt)  =  H'(R)e±ia)t 

then  the  probability  of  the  system  being  in  state  m  at 
time  t  (after  starting  in  state  k)  for  a  dipole  trans- 
ition is 

(7)     |Cm(t)|2  -  41^  |2sin[(Ek-  Em  ±ha>)  t/2h]  /  (Ek-Em±hto)  : 
where  Em,  Ek  are  the  energies  of  the  m  and  k  states. 
The  case  where  Em=  Ek±hoj  is  called  resonance  absorption 
or  emission  (depending  on  the  sign) ,  and  for  this  case 
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the  transition  probability  depends  linearly  on  time   for 

short  times.   (Recall  that  this  is  only  a  perturbation 

approach  and  the  interaction  must  still  be  small  or  occur 

for  a  short  time-  )   The  matrix  elements  coupling  the 

states  may  be  written 

H^  =  <  m|H-(R)|k  >  =  <Ym|ir  (R)  |Yk> 

where  the  wave  functions  are  those  calculated  from  the 

unperturbed  Hamiltonian.    A  first  order  calculation 

results  in 

Hrar  =  ^y  (mk) 
where  y (mk)  is  the  transition  dipole  moment 

y  (mk)  =  <  m|  y  !  k  > 

and  y  is  the  permanent  dipole  of  the  molecule. 

The  probability  of  a  level  changing  from  state  k  to 

m  in  a  radiation  field  with  a  finite  spectrum  rather  than 

simply  a  monochromatic  field  may  be  found  by  integrating 

(7)  over  the  frequency  domain  and  averaging  over  all 

spatial  directions.   The  result  may  be  written  in  terms 

of  the  Einstein  stimulated  emission  coefficient 

Pk-wa  =  UvBmk 
where  U   is  the  radiation  density  at  frequency  corres- 
ponding to  energy  difference  of  levels  m  and  k,  and  Bmk 
is  the  Einstein  stimulated  emission  coefficient 

Bkm  =  8u2/3h2 |y (mk) | 2 
Transitions  which  are  forbidden;  those  for  which  the 
transition  dipole  (or  dipole  moment)  are  identically  zero 
are  found  by  symmmetry  operations  and  it  is  usually 


34 

relatively  easy  to  find  which  states  are  not  infrared 

active.  33 

The  mechanism  by  which  laser  action  initiates  and 

subsequently  amplifies  is  most  easily  explained  with  use 

of  the  Einstein  coefficients.   In  the  early  part  of  this 

century  it  was  shown  by  A.  Einstein  that  a  resonant  photon 

travelling  in  one  direction  could  react  with  an  excited 

molecule  and  produce  two  coherent  photons  traveling  in 

phase,  plus  an  unexcited  molecule.   Once  laser  action  is 

started,  the  laser  medium  acts  like  an  amplifier.   The 

probability  of  an  excited  system  undergoing  such  a 

transition  in  a  radiation  field  of  energy  density  Uv  has 

already  been  given  as 

^2->i  =  UVB12 
and  similarly  for  absorption.   This  does  not, however, 

explain  how  any  radiation  gets  started  in  the  cavity  since 

initially  for  zero  radiation  there  is  zero  probability  of 

initiation.   Actually , atoms  or  molecules  in  excited 

states  remain  so  only  for  limited  periods  of  time  due  to 

the  Heisenberg  uncertainty  principle.   The  uncertainty 

may  be  expressed 

AE  At  -    h      or 

Av  At  a  1 
so  that  for  an  ensemble  of  molecules  of  lineshape  with 
half  width  Av  the  excited  state  can  be  expected  to  decay 
spontaneously  within  At.   Therefore  we  must  add  a  term  to 
our  probability  for  a  molecule  to  go  from  state  2  to 
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state  1:     P2^  =  UVB21  +  A21  . 

A21  is  the  Einstein  spontaneous  emission  coefficient  (At"1) 
Thus  we  see  that  it  is  through  the  natural  decay  that 
noise  is  able  to  start  the  laser  amplifier. 

Vibrational  or  rotational  energy  can  be  transferred 
collisionally  as  well  as  through  radiative  transition 
(for  infrared  we  do  not  consider  electronic  transitions)  . 
Collisional  transitions  occur  when  one  molecule  collides 
with  a  second  and  transfers  some,  or  all, of  its  energy  to 
the  different  degrees  of  freedom  of  the  second  molecule. 
This  energy  may  go  into  rotational,  vibrational,  or 
translational energy  of  the  second  molecule.   The  time  a 
molecule  remains  in  an  excited  state  before  the  energy 
redistributes  back  to  equilibrium  conditions  is  called 
the  relaxation  time.   In  optical  pumping,  energy  transfer 
can  proceed  three  ways;  transfer  from  one  excited 
vibrational  state  to  another  (V-V) , transfer  from  an 
excited  vibrational  state  to  higher  rotational  energy 
(V-R) , and  transfer  from  an  excited  vibrational  state  to 
the  translational  motion  of  the  molecule  as  a  whole  (V-T) . 
It  is  commonly  believed  that  tv_v<  tv_r<  tv_t  where  x 
is  the  relaxation  time  for  that  particular  process  and 
we  assume  that  transfer  within  one  vibrational  mode  is 
faster  than  between  two  different  modes.20'23'35 

Theoretical  accounts  of  collision  processes  have  only 
moderate  success  in  predicting  probability  of  energy 
transfer.   The  reader  is  referred  to  an  excellent  review 
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of  the  theoretical  treatment  of  collisional  energy  transfer 
by  Bailey  and  Cruickshank.36  The  most  popular  theory  to 
compare  with  experimental  results  is  a  quantum  mechanical 
V-T  theory  developed  by  Schwartz,  Slawsky,  and  Herzfeld 
(SSH  theory).37  SSH  theory,  as  well  as  most  others,  does 
not  succeed  in  providing  absolute  probabilities.   However 
it  does  predict  that  resonant  processes  (those  for  which 
the  second  of  the  colliding  pair  has  an  energy  level 
resonant  with  the  excited  level  of  the  first)  have  the 
largest  relative  probability  , and   that  pairs  with  a 
smaller  reduced  mass  will  have  higher  probability  of 
energy  transfer  than  others. 

Although  there  are  no  rigorous  selection  rules  for 
molecular  collisions,  there  appear  to  be  some  collisional 
transfer  processes  which  occur  with  much  higher  proba- 
bility than  others.   This  is  based  on  the  fact  that 
symmetric  states  do  not  combine  with  antisymmetric  for 
any  kinds  of  transitions  including  collisions.   It  is  for 
this  reason  that  the  two  modifications  of  symmetric  top 
molecules  such  as  NH3,  CH3C1,  CH3F  (e,a)  are  transferred 
into  each  other  only  extremely  slowly,  just  as  are  ortho- 
and  para-  hydrogen.  38 

Recently  Oka  has  given  some  approximate  selection 
rules  observed  by  microwave  studies  of  NH 3  and  the 
qualitative  interpretation  of  them. 39  To  summarize  these 
he  finds 


37 
(i)   Collision-induced  transitions  with  dipole  selection 
rules  (AJ  =  0,  ±1,  parity  +  ++-)    are  "preferred", 
(ii)   The  AJ  =  0  dipole-type  transitions  (same  selection 
rules  as  dipole  transitions)  have  much  greater  proba- 
bility than  the  AJ  =  ±1  dipole-type  transitions  for  levels 
with  J-K  but  they  have  probabilities  of  equal  order  of 
magnitude  for  levels  with  J>>K. 

(iii)   The  AJ>1  transitions  have  much  smaller  probabilities 
than  the  A J  =  ±1  transitions. 

(iv)   The  AK  /  0  transitions  have  much  smaller  proba- 
bilities than  the  AK  =  0  transitions. 

(v)  It  is  suggested  that  AJ  =  ±1  quadrupole- type  trans- 
itions (parity  +  -*-*•  +  )  are  the  same  order  of  magnitude  as 
those  of  the  corresponding  dipole  transitions. 

These  "selection"  rules  are  very  important  since 
most  experiments  measure  the  rate  of  transfer  of 
vibrational  energy  between  different  modes  of  a  molecule 
or  different  molecules,  although  they  don't  give  information 
as  to  the  particular  rotational  levels  involved. 

We  have  seen  in  this  chapter  that  there  are  two  basic 
energy  transfer  mechanisms:  radiative  and  collisional. 
Both  are  important  to  lasers,  sometimes  beneficial  some- 
times detrimental.   Radiative  transitions  are  obviously 
most  important  when  optically  pumping  a  gas  with  a  laser 
such  as  C02;  however  if  laser  action  is  expected  from 
any  level  other  than  that  pumped,  one  must  rely  on 
collisions  to  transfer  energy  selectively  into  that  level. 
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Transitions  due  to  the  natural  lifetimes  of  the  state 
must  initiate  stimulated  emission  before  collisions 
populate  the  lower  state,  thus  depleting  the  population 
inversion.   Both  spontaneous  and  stimulated  coefficients 
must  be  known  as  well  as  the  collisional  rates  and  their 
selection  rules  if  one  is  to  successfully  predict  how 
a  collection  of  molecules  will  react  after  optical  pumping. 
The  regeneration  of  laser  radiation  at  a  different 
frequency  will  be  a  function  of  all  three  processes. 


CHAPTER  III  -  METHYL  FLUORIDE 

The  success  of  producing  an  optically  pumped  laser  is 
largely  dependent  on  the  gas  to  be  used  as  a  lasing 
medium.   It  must  absorb  C02  laser  radiation  strongly 
and  be  able  to  achieve  a  population  inversion  between  two 
levels  that  would  result  in  16  micron  emission.   When 
modeling  a  process  such  as  this,  a  good  deal  of  informa- 
tion is  required  to  accurately  predict  the  results. 

The  first  area  that  must  be  understood  is  the  absorp- 
tion process.   It  must  be  known  how  strongly  and  on  what 
C02  line  the  gas  absorbs.   A  high  resolution  spectrum  is 
helpful  to  see  to  what  degree  the  gas  and  pumping  line- 
shapes  overlap.   Normally,  one  is  not  given  this  informa- 
tion, rather  an  experimental  absorption  coefficient  is 
given  for  a  particular  C02  line.   The  overlap  integral 
may  still  be  estimated  if  the  Einstein  coefficients  can 
be  found  from  vibrational  band  intensity  measurements 
which  have  been  carried  out  for  many  molecules.   In 
general  the  absorption  coefficient  will  depend  on  the 
particular  rotational  levels  involved,  line  broadening, 
degree  of  resonance,  and  pressure.   One  advantage  of 
modeling  a  process  such  as  this  is  that  these  parameters 
may  be  changed  and  the  absorption  coefficient  still  known 
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if  a  measurement  has  been  carried  out  for  any  one  partic- 
ular   set  of  parameters. 

The  second  area  of  importance  is  the  collision 
kinetics  of  the  molecule.   This  is  important  to  the 
problem  of  rotational  bottlenecking.   Often  during  laser 
operation,  a  particular  rotational  state  will  be  filled 
faster  than  the  collisional  rotational  relation  can  relax 
it.   When  this  occurs,  the  population  inversion  cannot 
be  maintained  and  lasing  stops.   The  bottleneck  effect 
is  so  termed  because  it  is  normally  the  process  that 
limits  laser  efficiency.   Within  this  area  also, is  the 
V-V  energy  transfer  process.   If  laser  action  is  expected 
from  any  level  other  than  the  upper  state  of  the  pumping 
process,  energy  will  have  to  be  transferred  via  V-V 
collisions.   The  time  scale  and  efficiency  of  these 
collisions  will  be  very  important  to  the  success  of  the 
laser. 

The  gas  to  be  used  in  generation  of  a  16  micron  laser 
by  optically  pumping  with  C02  must  strongly  absorb  C02 
radiation,  and  must  have  two  levels  separated  by  approx- 
imately 624  cmT1   The  lower  level  must  be  far  enough 
above  ground,  that  there  is  no  appreciable  thermal 
population.   Although  several  molecules  fulfill  the 
above  requirements,  only  one  has  had  enough  experimental 
work  done  on  it  to  make  it  attractive.   The  molecule, 
methyl  fluoride  (CH3F) ,had  been  used  as  a  saturable 
absorber  in  CO?  mode-locking  studies  and  far  infrared 
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laser  generation ,  30" 32  a  source  of  photon  echoes  in  phase 
coherence  studies , k0   and  recently  the  collisional  kinetics 
have  been  studied  by  following  laser  induced  f luoresence. 25 

Figure  15  displays  the  energy  level  diagram  for  methyl 
fluoride.   The  ground  state  to  V3=  1  transition  is  at 
approximately  9.55  microns  which  is  within  the  9.6  band 
of  C02."1  The  second  excited  state  of  V3  is  reported  at 
2081  cm"1  and  the  lower  pair  of  the  V25  doublet  hybrid 
band  is  reported  approximately  at  1460  cm-  . 

Methyl  fluoride  is  a  symmetric  top  molecule  with  C3v 

symmetry.38  The  rigid  rotor  term  values  are  given  by 

F(J,K)  =  AJ(J+1)  +  (A-B)K2 

where  J  is  the  total  angular  momentum  and  K  is  the 

component  of  J  on  the  internuclear  axis.   The  constants 

A  and  B  are  defined  as 

A  =  h2/2cla  ,   B  =  h2/2db 

where  I  is  the  moment  of  inertia  about  that  particular 

axis.   Selection  rules  yield  radiative  transitions  (rigid 

rotor,  harmonic  oscillator)  AJ=0 , ±1 , AK=0 , ±1  which  lead 

to  the  following  term  symbols,  labeled 

KJ(Ji     ,K,     ) 
lower   lower 

Qq(j,k)=v0 

Qp(J,K)=v0+2AJ 
Qr(J,K)=v0-2A(J+1) 
PQ(J,K)=V0-(A-B) (2K+1) 
PP(J,K)=v0+2AJ-(A-B) (2K+1) 
PR(J,K)=V0-2A(J+1)-(A-B) (2K+1) 
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Figure  15. Energy  levels  of  CH3F 
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RQ(J,K)=V0+(A-B) (2K+1) 

RP(J,K)=v0+2AJ+(A-B) (2K+1) 

RR(J,K)=v0-2A(J+l)  +  (A-B)  (2K+1) 
Actually  the  rotational  constants  change  slightly  from 
one  vibrational  level  to  the  next  and   centrifugal 
distortion  terms  may  be  included  so  that  the  spectrum  is 
slightly  more  complicated  than  is  indicated  here. 

Thus  it  may  be  seen  that  the  spectrum  of  methyl 
fluoride  is  quite  complex.   Parallel  bands,  those  for 
which  the  dipole  changes  parallel  to  the  applied  field 
have  selection  rules  AK=0,  while  perpendicular  bands 
have  AK=±1.   It  is  customary  to  make  the  following 
vibrational  assignments  for  methyl  fluoride."3  Actually 
the  V2  and  V5  vibrations  must  be  considered  degenerate 
thus  creating  a  V25  and  V52  hybrid  pair.   This  band  is 
composed  of  both  a  perpendicular  and  a  parallel  band. 
The  transition  from  V3  =  2  to  V25  to  be  considered  should 
therefore  have  several  rotational  lines  on  which  lasing 
could  occur  around  16  microns. 


Vibration 

Frequency (cm-1 

Species 

Type 

VJ-CH3    s-stretch 

2930 

a 

parallel 

V2-CH3    s-deform 

1464 

e 

parallel 

V3-  CF    stretch 

1049 

e 

parallel 

Vk-  CH3    d-stretch 

3006 

e 

perpendicular 

V5-  CH3    d-deform 

1467 

e 

perpendicular 

V ,-CH,    rock 

1182 

e 

perpendicular 

4  4 
Methyl  fluoride  has  been  used  as  an  absorber  of 
9.6  micron  C02  radiation  and  absorption  on  the  P(20) 
and  P(32)  CO,  lines  has   been  reported.   The  P(20) 
absorption  is  believed  to  be  the  QQ(12,1)  and  Qq(12,2) 
meythl  fluoride  transition  while  the  P(32)  absorption  is  due 
to  another  isotopic  species  of  methyl  fluoride.30"32  The 
lack  of  absorption  of  other  C02  wavelength  is  somewhat 
suprising  since  a  medium  resolution  spectrum  exhibits  a 
large  Q  branch  centered  around  the  P(18),  9.6  C02  band.1*1* 
Closer  inspection  reveals  the  K  splitting  of  this  branch 
may  account  for  the  possibility  that  the  extremely 
narrow  laser  line  can  fit  between  two  methyl  fluoride 
transitions.   The  laser  used  for  this  experiment  was  a 
Q-switched  mode-stabilize  C02  and  the  possibility  of 
hole  burning  exists.   The  effect  of  a  TEA  laser  pulse  or 
a  mode  stabilized  TEA  laser  can  be  modeled   as  described 
in  Chapter  I.   This  can  be  written  into  any  simulation 
scheme  and  results  checked  to  observe  the  effect  of 
different  forms  of  pumping  radiation. 

Recently  the  collisional  kinetics  and  energy  transfer 
have  been  studied  by  monitoring  fluorescence  rise  times 
from  various  levels  after  pumping  the  V3  =  l  with  a  C02 
laser .  25' 35' lt5~'+8  This  technique  has  been  fruitful  in 
observing  collisional  energy  transfer  from  a  particular 
vibrational  level  to  a  diluent  gas  also.1*9  Figure  16 
indicates  the  relative  speed  of  V-V  transfer  after 
pumping  the  V,=l  level  in  methyl  fluoride.   The  results 


45 


(1)  Excitation  of  V3  by  9.6  P(20)  CO2        _1 

(2)  2CH3F(V3)^r  CH3FIO)  +  CH3F(2V3)  +  10  cm  _1 
CH3F(2V3)  +  CH3F(V3)-^CH3F(0)  +  CH3F(3V3)  +  20  cm 

10  collisions  _j 

(3)  CH3F(3V3)  +  CH3F(0)^-  CH3F  (0)  +  CH3F(Vi)  +  120  cm_  ± 
CH3F(3V3)  +  CH3F(0)^r^-  CH3F  (0)  +  CH3F(V4)  +  100  cm 

70  collisions 

(4)  CH3F(Vi)  +  CH3F(0)^-  CH3F  (0)  +CH3F(2V2) 
CH3F(V!)  +  CH3F(0)^=r  CH3F(0)  +  CH3F(2V5) 
CH3F(V4)  +  CH3F(0)^=^  CH3F(0)  +CH3F(2V2) 
CH3F(V4)   +  CH3F(0)^-  CH3F(0)  +  CH3F  (2V5) 


10  collisions 


■10    cm 


•1 


(5)  CH3F(2V2)     +    CH3F(0)^^  2CH3F(V2;     - 
CH3F(2V5)    +    CH3F(0)  — 2CH3F(V5)-10    cm 

50   collisions 

(6)  CH   F(V3)       +    CH3F  (0)  ^^  CH3F(0)     +   CH3F(V)-133    cm 

40  collisions 

(7)  CH3F(V-T/R)   15,000  collisions 


Figure  16.  V~V,V-T/R  rates  for  methyl  fluoride. 
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were  obtained  by  monitoring  the  fluoresence  risetime  in 
an  experiment  similar  to  one  described  in  Chapter  I. 
The  processes  are  step-wise  in  that  they  follow  in  the 
order  indicated.   Comparison  with  SSH  breathing  sphere 
amplitudes  reveal   the  relative  collision  rates  are  in 
good  agreement  although  as  is  usually  true,  absolute 
probabilities  are  not."6'"7 

The  very  long  V-T  time  is  due  to  the  fact  that  the 
lowest  level  of  methyl  fluoride  is  still  over  1000  cm"1 
above  ground.   By  using  a  double  resonance  experimental 
setup,  where  a  particular  rotational  level  is  populated 
into  the  upper  vibrational  state,  and  microwave  absorption 
monitored  as  a  function  of  time   for  AJ=1,  AK=0  in  C13  H3F 
the  rotational  relaxation  was  found  to  be 

T  =  (10.5±6)ysec/mtorr.  which  is  long  as  compared  with  the 
V-V  rates.  50 

As  can  be  seen  from  the  preceding   discussion,  a 
great  deal  is  known  about  methyl  fluoride.   Figure  17 
indicates  the  match  up  between  known  C02  laser  lines  and 
the  theoretical  methyl  fluoride  spectrum  around  the 
reported  P(20)  absorption.^'51  Other  C02  lines  have  close 
coincidence,   although  we  will  primarily  discuss  the  P(20). 
This  figure  clearly  indicates  that  even  high  resolution 
is  not  enough  to  guarantee  good  absorption  of  C02 • 
Due   to  the  number  and  narrowness  of  methyl  fluoride 
absorption  lines,  one  might  expect  little  absorption  of 
a  CO,  line  that  chanced  to  fall  between  two  methyl 
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j     k ^Q  term  value 

"l2      "   0  1046.837*"    Reported  C02  line  1046.85 

12                        1  1046. 82l\             P(20)  absorption 

12          2  1046.826^      reported. 

12          3  1046.824 

12          4  1046.821 

12          5  1046.817 

12          6  1046.813 

12          7  1046.808 

12          8  1046.802 

12          9  1046.796 

12          10  1046.788 

12          11  1046.780 

12          12  1046.771 

Half  width  CH3Fr  2.21  x 

Half  width  C02;  2-32^x 
Cavity  width  CO2 ;  10~Jcm 


Figure  17.  Theoretical  spectrum  of  CH3F 
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fluoride  lines.   The  question  of  how  pressure  broadening 
of  both  C0?  and  methyl  fluoride  affects  absorption  has 
not  been  answered.   It  is  probably  true  that  the   Q(12) 
transition  is  involved  instead  of  the  next  closest  QP(1) 
since  the  thermal  population  peak   is  around  J  =  11.  ' 
Lasers  which  operate  on  a  single  axial  mode,  tunable  over 
their  doppler  width  have  recently  been  made  commercially 
available  and  will  no  doubt  play  an  important  role  in 
studying  and  maximizing  absorption  in  molecules  such  as 
methyl  fluoride. 

Of  course   the  real  interest  of  this  study  is  the 
generation  of  a  16  micron  laser.   Figure  18  shows  the 
approximate  expected  spectrum  for  V3=  2  to  V25  transition. 
The  rotational  constants  are  not  accurately  known  for 
either  level  so  the  ground  state  to  V3=l  rotational 
constants  were  used.32'43  These  may  be  in  error  by  several  wave 
numbers;  however  it  serves  to  give  an  idea  of  the  relative 
differences  among  transitions. 

Although  this  paper  is  concerned  with  the  computer 
simulation  of  laser  construction,  Figure  19  illustrates 
a  typical  laboratory  construction  of  an  optically  pumped 
laser.   The  KBr  prism  acts  as  a  dispersing  element  thus 
enabling  reflection  of  the  C02  and  16  micron  radiation 
to  physically  take  place  in  different  regions.   The 
grating  is  needed  to  operate  on  the  9.6  band  since  gain 
is  normally  higher  in  the  10.6  region.   Mirrors  M,  and  M2 
should  be  gold  or  dielectric  coated  to  reflect  at  16  microns 


49 


J'=12,  K'=2 


J'=12,  K'=l 


F(J,K) 


F(J,K)    J"    K" 


^Q 

619.23 

1.2 

2 

619.23 

12 

1 

3R 

641.07 

11 

2 

599.07 

13 

1 

qp 

599.07 

13 

2 

641.07 

11 

1 

rQ 

631.57 

12 

2 

622.78 

12 

0 

?Q 

596.37 

12 

3 

605.17 

12 

2 

rR 

651.73 

11 

1 

642.93 

11 

0 

Pr 

616.52 

11 

3 

625.32 

11 

2 

rp 

607.73 

13 

1 

600.93 

13 

0 

p? 

674.53 

13 

3 

583.33 

13 

2 

Figure  18.  2V3-V25  approximate  spectrum. 


5  0 


jj 

■■;• 
S 

m 

o 

c 
o 

■H 
+J 

rJ 

+J 

LQ 
C 

o 

u 


=  e 

OD  •  — 


51 

If  a  line  other  than  that  at  16  microns  has  a  lower 
threshold,  it  will  lase  first.   In  this  event,  mirror  M2 
will  have  to  be  replaced  with  a  grating  so  that  the 
cavity  becomes  tuned  to  16  microns. 

It  can  be  seen  in  figure  19  that  there  are  many 
different  external  conditions  under  which  lasing  may  be 
attempted.   These  include   gain  length,  pressure,  mirror 
reflectivities,  pump  power  and  duration,  and  type  of  C02 
pumping.   The  purpose  of  this  simulation  is  to  check  the 
effect  of  different  parameters  and  narrow  the  conditions 
under  which  one  should  expect  laser  action.   This  will 
eliminate  much  trial  and  error  experimental  work  and 
should  help  in  the  planning  and  construction  of  the  laser. 

A  second  area  where  computer  simulation  is  expedient 
is  the  determination  of  how  accurately  a  particular 
molecular  constant  is  needed  to  be  known  for  accurate 
prediction  of  laser  generation.   This  will  help  plan 
additional  experiments  which  need  to  be  performed  prior 
to  laser  construction.   In  the  case  of  methyl  fluoride 
there  are  two  areas  of  importance  that  may  need  further 
study.   The  first  is  the  absolute  frequency  of  the 
V3=  2  to  V25  transition.   This  will  have  to  be  measured 
if  it  is  deemed  important.   The  second  area  of  importance 
is  the  life-time  of  the  V3=  2  state  in  methyl  fluoride. 
This  is  important  because  a  long  lifetime  (determined  by 
the  Einstein  spontaneous  emission  coefficient)  will  allow 
the  level  to  collisionally  populate,  thus  depleting  the 
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population  inversion.   The  spontaneous  emission  acts  as 
an  initiator  of  laser  action  and  must  initiate  rapidly 
enough  to  allow  the  energy  to  transfer  radiatively  from 
V  =  2  to  V25  =  1  rather  than  collisionally .   Thus  the 
competition  between  collisional  and  radiative  lifetimes 
may  be  extremely  important.   An  important  aspect  of 
modeling  as  we  have  described,  is  that  the  collision 
rates  are  known.   This  should  enable  us  to  predict  for 
what  range  of  values  of  the  Einstein  coefficient  to  expect 
lasing.   It  is  conceivable  that  a  rough  estimate  within 
a  factor  of  100  is  all  that  is  needed  to  assure  lasing, 
or  an  accurate  measurement  may  have  to  be  made  if  the 
results  are  very  sensitive  to  the  actual  numerical  value. 
We  believe  that  the  Einstein  coefficient  for  the  V3=  2 
state  may  be  estimated  fairly  accurately  from  the  Einstein 
coefficients  for  the  V3=  1  level.   This   coefficient 
may  be  determined  by  intergrated  intensity  measurements 
using  the  relation 

r(cm2  mole"1)  =  2.505  x  102b{2 

52 

substituting  the  measured  value  for  T   yields. 

r  =  9055.6  ±10% 

B*2  =  6.8  x  10 23  an3  J-1  sec-2 
It  is  expected  that  the  Einstein  coefficient  for  the  V3=  2 
to  V25  =  1  will  be  substantially  smaller  than  B12  since  the 
transition  is  forbidden  by  harmonic  oscillator  selection 
rules  (AV  f    2)  .   Thus  by  estimating  the  coefficient  for 
this  transition  we  may  determine  roughly  to  what  accuracy 
it  must  be  known.   See  Appendix  B. 


CHAPTER  IV  MODEL 

The  model  which  we  employ  is  a  standard  kinetic 
rate  equation    used  by  many  authors.   The  two  types  of 
terms  are  those  for  radiative52  and  collisional 
transitions.    Figure  20  shows  the  interaction  of 
different  levels  in  this  model.   This  figure  should  be 
compared  with  the  energy  level  diagram  for  methyl 
fluoride  in  figure  14.   We  assume  that  collisionally 
AK  =  0,  AJ  =  0  selection  rules  predominate  so  that 
relatively  few  rotation  levels  need  to  be  represented. 
The  V25  state  is  represented  by  two  rotational  levels, 
one  for  collisional  population  and  one  for  radiative 
population  since  it  is  desired  to  see  what  effect  the 
different  symmetry  combinations  have  on  lasing. 
Rotational  levels  of  methyl  fluoride  are  of  two  types, 
the  K  =  0,3, 6,... are  symmetric  (type  a)  while  the 
K  =  1,2 ,4,5, .. .levels  are  asymmetric  (type  e ) .   Levels 
with  overall  species  e  are  doubly  degenerate,  thus 
direct  product  tables  show  that  the  V,  level  has  rotational 
levels  a  and  e,  while  the  others  have  only  the  a 
rotational  level;  these  yield  the  overall  e  state 
(a  x  e  =  e)  while  the  e  rotational  level  gives  an 
a  +  a  +  e(exe  =  a  +  a  +  e).38  Molecules  of  methyl 
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Figure  20.  Schematic  representation  of  model. 

(  large  numbers  indicate  number  of  collisions 

and  arrows,  the  direction  for  that  rate.) 
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fluoride  in  the  ground  state  are  either  in  rotational 
type  a  or  e  so  that  even  collisionally ,  little  mixing  is 
expected.   The  possibility  of  changing  rotational  levels 
by  AK  =  0,  +  l  exists,  since  V2  is  actually  a  hybrid 
(parallel  and  perpendicular  allowed  transitions)  and 
thus  the  lower  laser  level  could  possibly  be  a  different 
symmetry  than  the  upper  level  which  is  filled  collisionally. 
The  equations  when  non  dimensionalized  in  time  with 
respect  to  collision  frequency,  number  density  with 
respect  to  density  of  the  ground  rotational  level,  and 
laser  output  with  respect  to  the  input  C02  power  assume 
the  form 


51  =  0¥  (N    -  g?v^   NL  )  +  PH< 
dt   y  u  2V3   g„,c   V25 


^  =  r2  (N0NW,  -  e1Q/kT  N23  )  +  r,  (N0N3V3  -  e20/kT  N2v,Nv3) 
+  r12  (N0NV5  -  e133/kTNo  Ny  s ,  +  ZU*(N0-  Ny3  ) 

3j|?  =  r2(e^kT  N23-  N0N2V3)  +  r,(N0N3V3-  e2°/kT  N2v3NV3) 
Z23  ^(N2V3  -  SL>^   nLc) 

V25 


|^  =    r,(e20'kT    N2v3NV3   -  N„   N3v3 )    +  r5(N0Nvltl-    e110/^    N0N3v3) 


dNv*i  =    r     (ello/kT    N   N         -   N   N         )    +    r,{NnN,„,q     -   NnN.rui   ) 

— jl  ■l5v  0       3V3  0      V  h  1  6V       0       2  V  2  5  0      V*  1     ' 
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lP      -    r6(N0Nv41-   N0N2V25)    +    r7(N«5     -    e10  ^  N0  N  2V25 ) 


dN 


dt 


+    ra(NA        "   e10/kT    NN,VJ 


r7(e10/kT     N0N2V25      -   N^5  )     +    r9  (NQ  NV52     -N0N$a5> 


+   rF   Z„  ¥ 


"    r10N0NV25  +  rnN0NV6    -    TNV25 
(N2V3-    |^L3     n£2S  )-    U(Nj25) 

"v25 


^        -    Z23^(N2V3-    |2L3     N^25)-TNt25     +rF[r7(e^T     NQN2v25 

y  2V3 

-   n^225)    +    r9(N0NV52  -   NX25  )-    r10(N0N^25  ) 
+    r„  N0  NV6  ]-    U(N^2.  ) 


dN^.        =    r8(e10/kT    Nq  N^v25      -   Nv252  )    +    r9(N0<25     -   NflNV52  ) 


"    r10NQNV52    +    ru  N0NVR     -    YNV52 


^l±        =   r12    (e113^1  N0NV3     -    N0NV5)     -    2rn    N0  N^ 
dt 

+    r10N0(NV25     +    NV52  )     -    XNV6 

N0    =    K    -     [NV3    +   N2V3     +    ^V3      +    N2V25     +   N^25 
+    NL         +    N  +    N„    ] 

V2  5  V5  2  V6 

where 

¥  =  16  micron  laser 

$  =  input  C02  laser 

N.=  number  density  of  particular  energy  level, 
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where 

9 

K  «£  N.(t=0) 

i=o 

and  we  define  the  following  constants 

Q   =   hv0N0B23  W 

v 

P   =   C    In   r1  r2 
2Lv 

Z        =     0,     V^    B12  V 

Vq  CV 

z23=  fi  ^i  b^w 

v  u       cv 

r.    -   k.  N°/v 
i  i 

(Collisional  reaction  rate  k..  is  defined  by 

dNi,  =  v.  N.N-  -  k-  N-Ni 
dt  J 

for  a  process  such  as 

CH3F(Ni)  +  CH3F(Nk)  ;-  CH3F(Ni)  +  CH3  F  (Nj  )  +  *v 

from  the  equilibrium  condition 

Kji  ..  [CH3F(Nj)CH3F(Nj)]  =  ^v/kT  ) 
Ki7    [CH3F(Ni)CH3F(Nk)  J 

T  =  r/v  rotational  relaxation 

RF  =  %  relaxation  between  two  V25  rotational  levels. 

If  levels  are  of  different  species  RF  =  0 . 
U,Y,X  =  collisional  depopulation  with  diluent  gas  for 
that  particular  level. 

Initial  conditions  on  the  vibrational  level  popula- 
tions are  given  by  Boltzmann statistics ;  while  the  initial 
condition  on  the  16  micron  laser  is  given  by  the 
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spontaneous  emission.   The  following  parameters  are  used 

in  the  above  definition: 

h  =  Plancks  constant  =  6.626  x  10~27erg  sec 

v'  =  frequency  of  16  micron  Laser 

v0  =  frequency  of  C02  Laser 

B12  =  Einstein  stimulated  emission  coefficient  for 

ground  to  V3  transition. 
B23  =  Einstein  stimulated  emission  coefficient  for  2V3 

to  V25  transition. 
A23  =  Einstein  spontaneous  emission  coefficient  for 

2V3  to  V25  transition. 
N°  =  Population  of  the  ground  state  rotational  level 

which  is  optically  pumped. 

v  =   collisional  frequency 

W  =  methyl  fluoride  lineshape  at  16  microns 

V  =   overlap  between  pumping  lineshape  and  absorbing 

lineshape. 
Q,   -      Max  power  per  unit  area  for  C02  pumping  laser 

c  =   velocity  of  light  in  vacuum  3  x  10 10  cm  sec 
r  r   =  mirror  reflectivities 
1  =   laser  length  (cm) 

The  effect  of  tuning  the  laser  around  the  methyl  fluoride 
line  is  contained  within  the  constant  V.   From  the  quoted 
absorption  coefficient  the  overlap  intergal  may  be 
evaluated  numerically.3^32  Using   equation  Al ,  (Appendix  A) 
we  see  this  gives  a  distance  between  the  C02  and 
absorbing  methyl  fluoride  line  of  approximately 
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1.61  x  10" 2  cm-1.   This  agrees  favorably  with  the 
predicted  distance  from  the  theoretical  spectrum  of 
figure  17.   The  laser  used  in  this  experiment  was  a  Q- 
switched,  mode-stablized  laser.   By  changing  the  frequency 
difference  between  C02  and  methyl  fluoride  lines,  the 
effect  of  tuning  the  transition  to  exact  resonance  may 
be  simulated.   Figure  21  shows  the  predicted  position 
for  a  common  Q-switched  laser  relative  to  methyl  fluoride 
(Appendix  A) .   This  situation  may  be  changed  however  by 
changing  (v-vc).   An  experimental  condition  such  as 
figure  22  may  be  simulated  by  bringing  the  C02  laser  line 
very  close  to  the  methyl  fluoride  transition,  or  moving 
it  across  the  methyl  fluoride  transition  in  time. 

The  Einstein  coefficients  for  the  V3 =  2  to  V25 =  1 
transition  may  be  varied  to  see  how  sensitive  laser 
operation  is  to  them.   As  a  first  guess  (see  Appendix  B) 
we  choose  B23  -  4.58  x  1019  cm3  J_  J  sec" 2  .   The  collisional 
rates  are  given  in  figure  19  and  diluent  collision  rates 

j      1*5  ,1+g 

will  either  be  taken  from  experiment  or  assumed. 

The  rotational  relaxation  rate  of  methyl  fluoride, 
as  previously  stated,  has  been  measured  to  be  on  the  order 
of  103  collisions.   We  will  use  1000  collisions.   All 
calculations  are  carried  out  at  300  degrees  Kelvin.   Only 
these  parameters  which  may  be  easily  met  in  the  laboratory 
will  be  considered   thus  eliminating  exotic  experimental 
set-ups  to  check  the  computer  predictions. 


CO2  doppler 
Profile 


axial 
mode 


CH,F  J=12 


AG 


01  2      3 


Figure   21.    C02    and  methyl    fluoride    lineshapes. 
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Pressure   broadened 
C02 


Figure    22.    Broadened  C02    in   resonance  with 
methyl    fluoride. 


CHAPTER  V  RESULTS 

Numerical  solutions  of  the  equations  described  in 
the  last  chapter  were  carried  out  at  the  University  of 
Florida  on  an  I.B.M.  360-75  computer.   The  numerical 
technique  employed  was  Hamming's  modification  of  the  Milne 
predictor  -corrector  method.   A  fourth  order  Runge-Kutta 
method  is  used  to  generate  the  first  time-increment 
solution  since  the  predictor  -corrector  method  is  not 
self -starting.   Hamming's  method  is  a  stable  fourth 
order  integration  procedure  which  has  the  advantage  of 
a  variable  step  size.   This  saves  computing  time  without 
sacrificing  accuracy. 

The  accuracy  for  all  the  results  was  kept  between  one 
part  in  10 3  to  one  part  in  104.   This  is  better  than 
actually  required;  however  it  facilitates  faster  overall 
integration  since  solutions  were  not  allowed  to  start 
diverging  at  any  point.   The  time  step  was  between  .5 
and  .01  measured  in  units  of  time  where  one  unit  was  the 
time  for  one  collision.   The  predictor  -corrector  method 
was  able  to  bisect  the  time  step  up  to  10  times  if 
required  to  obtain  the  specified  accuracy. 

The  first  question  which  needed  to  be  answered  was 
what  parameters  to  start  the  model  with  since  there  are 
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B12  =  4.58  x  1019  cm3J_1  sec"2 
B12  =  4.58  x  1020  cm3J_1  sec-2 
Bia  =  4.58  x  101S  cm3J_1  sec"2 
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a  large  number  for  which  threshold  and  subsequent  lasing 
would  not  occur.   This  depends  on  the  Einstein  coefficient 
as  was  discussed  in  Chapter  IV.   Figures  23  and  24  show 
the  reflectivity  needed  at  a  particular  pressure  to 
satisfy  the  Schawlow-Townes  condition.   The  population 
inversion  may  be  estimated  from  the  thermal  population 
of  the  lower  laser  level  since  little  additional 
population  is  expected  during  the  first  several  collisions 
The  three  curves  are  drawn  for  the  Einstein  coefficients: 
(L  =  100cm) 

A, 
B. 
C; 

All  numerical  work  was  performed  for  case  A  although  if 
the  actual  B12  is  anywhere  in  the  above  range,  the 
equations  will  take  the  same  form  by  picking  a  point  on 
the  appropriate  curve.   The  laser  output  however   will 
decrease  as  the  reflectivity  increases  so  that  for 
reflectivity  above  .99  little  output  will  be  expected. 
It  can  be  seen  that  case  C  is  more  or  less  a  limiting 
case  for  practical  application. 

Unless  otherwise  stated  the  conditions  and  parameters 
were  as  follows: 

Pressure  =  10  torr 
Pump  power  =  10  watts/cm2 
Pump  duration  =  CW 
Reflectivities  =  100%,  98% 


6^ 


10 
Pressure  (torr) 


100 


Figure  23.  Schalow-Townes  condition  for  inversion  of 
5  x  10-3  (  L=  100cm  ) . 
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.96  ■ 


Pressure  (torr) 


Figure    24.    Schalow-Townes   condition    for    inversion   of 
1  x    10"  2,     (    L=    100  cm    )  . 
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Laser  length  =  200  cm 

Pump  J,K;  Jx=  12  Kl =    2 

Laser  transitions  (lower)  J,K;  Jn=  12,  Kn  =  2 

Distance  between  C02  and  CH3F  linecenters 

=  1.62  x  10" 2  cm"1 

Rotational  relaxation  =  10 3  collisions 

Diluent  relaxation  =  none. 
Cases  of  departure  from  these  conditions  will  be  dealt 
with  as  necessary.   Figure  25  shows  the  laser  pulse  as 
it  first  starts  to  develop  for  250  nsec  after  being  pumped 
by  a  kilowatt  C02  laser  pulse  for  50  nsec.   Notice  that 
the  laser  does  not  immediately  start  to  amplify  after 
threshold  is  reached,  but  waits  several  hundred  nanoseconds 
before  the  pulse  starts  to  curve  upwards  again.   This 
is  due  to  the  relatively  small  Einstein  coefficient 
C~1019  J"1  cm3  sec"2)  for  the  2V3  to  V?5  transition.   In 
studying  the  populations  of  each  level  in  the  model  during 
the  first  few  collisions  several  things  were  apparent. 
The  relatively  strong  absorption  coefficient  (.018cm  l torr 
for  P(20))  caused  the  ground  to  V3  transition  to  saturate 
within  a  few  collisions.   This  strongly  limited  the 
amount  of  energy  in  the  C02  pulse  which  was  utilized, 
since  after  saturation  nearly  all  the  laser  pulse 
propagates  through  the  laser  medium  as  a  bleaching  wave. 
Increasing  the  power  of  the  pump  C02  laser  for  this  type 
of  situation  was  found  to  have  little  effect  on  output 
or  populations   as  might  be  expected. 
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2  4 

collisions 

Figure  25.  CH3F  pulse  following  threshold. 
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In  the  figures  presented  for  energy  level  populations, 
the  following  symbols  are  used. 

SYMBOL  ENERGY  LEVEL 

2  V3 

3  2V3 

4  3V3 

5  Vw 

6  2V25 

7  vf5  collisionally  filled  level 


25 

9  V52 


V'L  laser  level 


0  V6 

Figures  26a  through  261  show   the  results  of  pumping 
methyl  fluoride  with  a  megawatt  lOOnsec  C02  laser.   The 
resulting  16  micron  laser  output  has  a  peak  intensity  of 
approximately  5  watts/cm'*  in  10" 6  sec.   Notice  that  only 
3.4  microseconds  elapse  between  pump  and  laser  output. 
The  maximum  population  inversion  occurs  in  only  20 
collisions  and  decays  nearly  linearly  for  the  next  150 
collisions.   The  population  of  V3=l  level  saturates 
immediately  and  decays  via  V-V  collisions.   The  populations 
of  the  V3=  1,2,3  levels  all  peak  with  a  phase  delay 
corresponding  to  the  V-V  equilibrium  rate.   Notice  that 
the  maximum  population  of  each  level  is  less  than  the 
level  previous  thus  indicating  the  collisional  equilibrium 
described  earlier.   This  is  not  true  for  the  V25  level 
which  is  also  populated  by  radiative  transitions.   It  is 


69 

interesting  to  note  that  although  the  V52  level  is  not 
populated  by  radiative  transitions,  it  is  always  in 
equlibrium  with  the  V25  level.   This  is  due  to  the  fact 
that  both  levels  form  a  hybrid  pair  and  are  different 
linear  combinations  of  the  same  two  levels.   These  levels 
are  degenerate  and  the  resulting  pair  are  in  Fermi 
resonance  and  are  thus  split  by  only  2cm-1.   In  contrast 
to  these  figures  consider  figures  27a  through  271.   This 
set  show  response  to  a  5  watt  CW  C02  laser.   Notice  that 
the  16  micron  laser  takes  nearly  2.5  times  longer  to 
form  its  pulse,  and  the  pulse  is  nearly  7  watts.   The 
long  phase  delay  is  understood  if  it  is  recalled  that 
lasing  occurs  from  the  V3=  2  level.   By  pumping  with  only 
5  watts  it  takes  nearly  100  collisions  to  saturate  the 
ground  to  V3  transition  and  thus  population  in  V3=  2 
does  not  peak  until  180  collisions  after  C02  pumping 
is  initiated.   The  V25  level  is  being  collisionally 
populated  during  this  time;  however  the  figure  clearly 
shows  that  the  rate  for  collisional  population  is  much 
smaller  than  that  for  radiative. 

This  explains  why  the  population  inversion  has  such 
a  slow  rise  time  and  rapid  decline,  as  compared  with  the 
previous  figure.   An  important  difference  to  note  between 
CW  and  short  pulsed  pumping  is  that  in  the  CW  case  the 
upper  vibrational  levels  are  collisionally  populated 
much  more  than  the  pulsed  case.   This  is  because  of  the 
long  phase  delay  of  the  output  laser  thus  requiring  levels 
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Figure  26a.  Parameters. 
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16  MICRON  LASER  OUTPUT ,HATT/CM2 
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Figure  26b.  Methyl  fluoride  laser. 
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Figure  26c.  Population  inversion, 
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Figure  26d.  Population  of  V3 
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Figure  26e.  Population  of  2V3 
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Figure  26f.  Population  of  3V3 
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Figure  2  6g.  Population  of  V41. 
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Figure  26h.  Population  of  2V25. 
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Figure  26i.  Population  of  V^5 
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Figure  2 6 j .  Population  of  V '   , 
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Figure  26k.   Population  of  V52 
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Figure  2  61.  Population  of  Vg, 
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Figure  27a.  Parameters, 
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Figure  27b.  Methyl  fluoride  laser. 
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Figure  27c.  Population  inversion, 


pmfPGV  LEVEL  POPULATIONS  ..-♦— + 


n 
) 

G 
R 

N 

0 

P 

n 
T 

L 

E 
V 
E 
L 


^'.500 


0.375 


O.250 


0.125 


222 


-?. 


".0 


222?2?222P 

T?2 


22 


2222 


2?222?2   I 

221 
I 

I 


0.0 


90.250 


1*0."500    270.750    3*1,000 


time 


Figure  27d.  Population  of  V3 


86 


ENERGY  LEVEL,  POPULATIONS  + 

O.250  +- -— * -♦  — -     ;        . 


K  ' 
Q 

) 

G 
R 

K 
0 

p 

n 

T 

L 
E 

v 

E 
L 


i 


0.188 


0.125 


0.063 


o.o 


I 

■531 
33   3 

,"       I 

T 
T 
I 

T 

T 

I 

T 
! 
I 
I 
I 
I 


333' 


1      i 

T    3333331 
33333333 


90.250 


l* 0.5 0  0 
time 


P7Q.7«;0    361.000 
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to  remain  in  equilibrium.   Thus  during  CW  operation  the 
lower  laser  level  (V25)  is  populated  much  higher  before 
lasing  thus  reducing  the  population  inversion.   The 
explanation  of  why  the  CW  case  has  higher  output  and 
lower  population  inversion  is  revealed  by  comparing  the 
two  figures  for  each  case  (figures  26c,  27c).   Although 
the  inversion  for  the  pulsed  input  reaches  a  higher 
value,  it  also  takes  nearly  three  times  as  long  to 
release  its  energy.   The  power  output  of  the  CW  pump  is 
then  understandably  higher  since  the  time  rate  of  change 
of  energy  released  is  higher. 

Of  the  22  different  parameters  sets  run  numerically, 
figures  26  and  27  are  representative  for  all.   As  might 
be  expected  there  was  no  magic  combination  of  parameters 
which  greatly  increased  laser  output,  although  all  22 
combinations  yielded  a  16  micron  laser  of  at  least 
several  watts/cm2.   The  figures  for  all  20  other  combina- 
tions will  not  be  shown  since  they  resemble  so  strongly 
those  of  26  and  27  that  little  additional  information 
would  be  gained.   (See  Appendix  C)  . 

Figure  28  indicates  the  efficiency  for  different 
types  of  carbon-dioxide  lasers.   The  parameter  I  is  a 
measure  of  energy  (I  =  .023mJ).   The  first  thing  one 
notices  about  this  figure  is  that  all  points  fall  more 
or  less  along  a  straight  line  which  peaks  around  log (I) =3 
and  sharply  drops  off  below  that.   This  represents  the 
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amount  of  energy  to  keep  the  ground  to  V3  transition 
saturated  while  the  other  vibrational  levels  are 
becoming  populated,  without  creating  a  bleaching  wave. 
Notice  that  in  all  cases,  where  pulse  time  lengthens, 
efficiency  increases.   As  you  demand  less  power  from  the 
system,  efficiency  is  able  to  get  quite  high.   The 
points  to  the  right  to  the  peak  have  much  higher  output 
power  than  those  to  the  left  ;  however  they  also  have 
much  lower  efficiency.   Pulse  times  for  all  points  are 
as  follows: 

tea;   100  nsec. 
q-switch;   1  microsec. 
CW;   CW. 
Quantum  efficiency  would  be  at  1.8  on  the  ordinate. 
Thus  it  is  seen  that  a  16  micron  laser  with  very  high 
efficiency  could  be  produced  if  one  were  content  to 
use  very  low  power  pump  C02  laser  and  get  low  power 
(several  watts  in  several  microseconds)  output.   This 
does  not  mean  however  that  high  radiation  density  cannot 
be  achieved  at  16  microns.   One  may  still  operate  a 
large  volume  methyl  fluoride  cavity  and  focus  the  beam 
to  increase  the  power  per  unit  area  at  the  focal  point. 


CHAPTER  VI  CONCLUSION 

The  feasibility  of  producing  a  16  micron  laser  from 
methyl  fluoride  has  been  demonstrated  in  this  paper. 
The  mechanism  for  energizing  the  methyl  fluoride  molecule 
is  through  optical  pumping  by  the  P(20),  9.6  band 
carbon-dioxide  laser  which  is  known  to  absorb  strongly 
into  methyl  fluoride.25'31'32  Numerical  work  was  done  on 
a  semi-classical  model  which  showed  this  technique  was 
indeed  possible.   The  success  of  this  type  of  approach 
is  due  to  the  large  amount  of  experimental  work  carried 
out  on  methyl  fluoride  previously,  thus  documenting 
absorption  coefficient,  absorbing  line ,  and  collisional 
rates .  25'  30'  n'  35'  w'  42'  4lf'  h5~  * 

The  methyl  fluoride  laser  described  herein  is 
characterized  by  high  reflectivity  and  pressure,  and 
long  gain  length.   This  is  expected,  since  the  amplifi- 
cation is  less  than  .51%/cm.  This  does  not  however 
preclude  the  methyl  fluoride  laser  from  having  substantial 
output  or  high  efficiency.   Output  of  from  5  to  20  watts/cm2 

in  several  microseconds  was  shown  to  be  obtainable 

under  several  conditions  and  configurations. 

Energy  storage  in  methyl  fluoride  is  small, 

(approximately  31.45  mJL_1torr_1)  so  that  methyl  fluoride 
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may  be  used  more  effectively  as  an  amplifier  of  16 
micron  radiation  than  an  oscillator.   This  would  be 
possible  since  small  amounts  of  16  micron  radiation  may 
be  generated  in  crystals.   The  system  envisioned  might 
be  similar  to  those  used  in  fusion  research.   A  small 
oscillator  would  feed  radiation  to  a  series  of  amplifiers 
which  would  provide  a  high  intensity  radiation  field. 


APPENDIX  A 


The  overlap  integral  between  the  C02  lineshape  and 
the  methyl  fluoride  lineshape  is  shown  below. 


CH3F 


C02 


AK, 


"o        "C 


Assuming  both  lines  are  homogeneously  broadened  the 
overlap  integral  takes  the  form 

v  =  fj$>   t(v-v0)2  +  (Avr,)2]"1  [(v-vc)2  +  (Av^)2]-^ 

let  to  =  v-v0   co1=vc-v0 

then  the  integral  becomes  V  =  AvcAym  j 

(2tt) 


where 


/ 


da) 


(a)2  +  A2)  ((co-oj1)2  +  B^) 


and 


A  =  ^m  ,   B  =  Av , 
2  2 


This  may  be  written  as  a  contour  integral  and  complete 
as  shown.  The  integral  has  two  simple  poles  at  iA  and 
iB  +  w1 . 
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From   the    figure 


/ 


dgj  =   i   +    r   =    2irilresidues 

(u/  +  A')  ((oj-u1)"  +  Bz) 


however  on  T    one  may  write 

u  =  reie 

and  the  integral  vanishes  as  r~ 3  as  r+~.   Thus  the 

contribution  along  r  =  0.   We  then  evaluate  the  residue 

at  each  pole  to  yield  the  result 

Al.    V  =  -L-[(^)2    +(^m)2+  (Ve-V^2]"1  (AvQ  +  AvJ 
2tt     2  2 

where  we  assume  that  the  frequency  distribution  is 

symmetric  with  respect  to  reflection  in  the  origin. 

This  is  assumed  because  photons  traveling  to  the  right 

through  the  cavity  see  the  frequency  shifted  up,  while 

those  traveling  to  the  left  see  it  shifted  down.   This 

is  only  strictly  true  for  inhomogenously  broadened  lines 

However  the  lineshape  between  the  two  is  similar  enough 

that  little  inaccuracy  is  introduced  by  assuming  a 

Lorentzian  lineshape. 


APPENDIX  B 

The  Einstein  coefficients  for  stimulated  and 
spontaneous  emission  from  the  V3  =  2  to  V2  or  V5  =  1 
levels  of  CH3F  could  be  obtained  from  the  measured 
integrated  molar  absorption  coefficients  from  infrared 
intensity  studies  of  the  "hot-band"  transition 
(V2  or  V5  =  1  to  V3  =  2) .   Unfortunately,  there  have 
not  been  any  such  studies  of  the  intensity  of  this 
transition.   Hence,  we  are  forced  to  estimate  the 
intensity  of  this  transition  from  some  theoretical 
treatment  of  "hot-band"  intensities.   The  most  under- 
standable treatment  of  this  problem  is  in  an  unpublished 
preprint  by  Yao  and  Overend.   They  find  that  the  matrix 
element  for  the  transition  dipole  for  a  difference  band  is 
<vs,vs.  ,vs..  |P'  |vs+l,vs.-l,vs..> 

v  +1   1/2   v  ,  V2      2P 
=  Pss'  -§ —       -§-     1   +   ^7  (g3(sss'_+  gjsss')} 


ss 


2P, 


s  ' 

+  p x    {g3(s's's)  +  gjs's's)} 

and  the  intensity  or  Einstein  coefficient  is  thus 

proportional  to  the  square  of  this  value.   Here  Pgs ,  = 

32P/30  30  ,  is  the  mixed  second  derivative  of  the  dipole 
s   s 
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moment  P  and  is  part  of  the  "electrical  anharmonicity" 
of  the  molecule,  s  and  s'  refer  to  the  vibrational 
modes  that  are  involved,  Ps  =  3P/3Qs-and  Pgl  =  3P/9QS- 
and  the  g   and  g k    functions  are  related  to  the  mechanical 
anharmonicities 

g3  (sss  '  )  =  03sus  ,q  1 

g4  (sss')  =  2w|ga 

g,  (sss')  =  -(ksss,/o,s)  (4o)|  -  <-)-1- 

Here  w   and  to  ,  are  the  harmonic  frequencies  involved, 

and  k    ,  is  the  mechanical  anharmonicity  potential  constant 

Since  we  do  not  know  enough  about  the  CH3F  molecule 
for  these  constants  (k    .  ,  P„  , ,  etc.)  to  be  available, 

SSS         •=»  23 

we  must  estimate  them.   We  note  that  the  particular 
difference  band  we  are  interested  in  is  the  transition 
from  v2(vs.)  or  v5  =  1  to  v3  =  2  (=vs  +  2).   Hence,  this 
theory  by  Yao  and  Overend  does  not  apply  strictly,  but 
we  may  expect  approximately  the  following: 

(1)  P   =0,  since  it  is  related  to  the 

i     intensity  of  an  overtone  transition 

forbidden  by  harmonic  selection  rules; 

(2)  P   ,  =  d3?/dq?3dQ2    =    0  since  the 
electrical  anharmonicity  is  small; 

(3)  ksss,/tos  =0.1,  so  that  g3+  g^'  0.1 

Hence 

<0,1,0|P|0,0,2>  s  0.1  (3P/9Q2)  . 
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We  conclude  that  a  first  guess  at  the  Einstein  coefficients 
for  this  transition (v3  =  2  to  v2  =  1)  is  that  they  are 
about   1   percent  of  the  Einstein  coefficients  for  the 
fundamental  (v2  =  1  to  0)  transition.   Similarly  the 
transition  from  v2  =  2  to  v5  =  1  is  about   1   percent 
of  the  value  from  the  v5  =  1  to  0  fundamental  transition. 
The  Einstein  coefficient  for  the  fundamentals  can  be 
calculated  from  the  measured  infrared  intensities  of  the 
fundamental  transitions,  summarized  for  example  by 
Russell,  Needham  and  Overend.52 
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APPENDIX  C 
PARAMETER  CHANGED 


OUTPUT  W/cm 
HALF  WIDTH 
(COLLISIONS) 


COMMENTS 


co2-103w/cm2 
4  collisions 


CO2-105W/cm2 
4  collisions 

CO2-106W/cm2 
4  collisions 

CO2-103W/cm2 
43  collisions 

CO2-105W/cm2 
43  collisions 

CO2- 10W/cm2 
CW(361  collisions) 

CO2-50W/cm2 
CW(361  collisions) 

CO2-100W/cm2 
CW(361  collisions) 

C02-lW/cm2 

CW(361  collisions) 


4.lW/cm2       1  collision 
90  collisions  =23  nsec. 
d>  10  torr 

4 .  6W/cm2 

80  collisions 

4 . 6W/cm2 

75  collisions 

1 6 . 4W/cm2      shor  te  r 

45  collisions  phase  delay 

16.5W/cm2 

45  collisions 

10 . 9W/cm2 

90  collisions 

15.2W/cm2      shorter 

55  collisions  phase  delay 

16.lW/cm2 

45  collisions 

no  output 


10 

C02-5W/cm2 

CW(361  collisions) 

7 . 6W/cm2 

130  collision 

longer 
tail 

11 

CO2-106W/cm2 
43  collisions 

16.5W/cm2 

45  collisions 

12 

L=100  cm. 

6 . 5W/cm2 

50  collisions 

compaire 
#12-22  with 

#6. 
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NUMBER 


PARAMETER  CHANGED 


OUTPUT 


COMMENTS 


13 
14 
15 
16 
17 
18 

19 

20 

21 

22 

23 


L=300  cm. 

r]_=1.0 
r2=.97 

r]_=1.0 
r2=.99 

CO2-CH3F  frequency 
difference=0 

J"  =  12,K"  =  1 
J"=12,K"=3 

pressure=25  torr 


diluent  relaxation 

of  V25-V52  in 
1600  collisions 

diluent  relaxation 
of  V,  in  1600 
collxsions 

diluent  relaxation 

of  V25-V52  in 
250  collisions 


13.4W/cm2 

110  collisions 

12.4W/cm2 

65  collisions 

7 . 6W/cm2      longer 
125  collision  tail 

16.8W/cm2 

65  collisions 

10 . 9W/cm2 

90  collisions 

16.5W/cm2      different 
75  collisions  species; 
K'=2,K"=3 

28.5W/cm2      1  collision 
200  collision  =9.2  nsec. 
(§)  25  torr 


diluent  relaxation 

of  V25-V52  in 

1600  collisions, 

Q-switch  in  100  collisions 


8 . lW/cm2 
8 . OW/cm2 
9 . 3W/cm2 
10.lW/cm2 


CW- lW/cm2 

CW-.52W/cm2 

CW-3.8W/cm2 

CW-3.8W/cra2 
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